use "C:\Users\AJUSTSA_G_VERDUGO\Documents\AJUSTSA_G_VERDUGO\Documents\ImmPanelRevis17\donnees\permanent_revis3y.dta", clear
keep if homme==1
gen first=(age>=25 & age<=50)
egen time=group(an)
egen nindiv=group(nninouv)
drop if missing(nindiv)
xtset nindiv time
tab first an
gen twice=(first==1 & !missing(F.time))
tab twice an if first, col

/* annual attrition rates */
cd "C:\Users\AJUSTSA_G_VERDUGO\Documents\AJUSTSA_G_VERDUGO\Documents\"
/* annual attrition rate */
use "ImmPanelRevis17\donnees\indiv_dispm.dta", clear
/* real wages */
joinby an using	 "C:\Users\AJUSTSA_G_VERDUGO\Documents\AJUSTSA_G_VERDUGO\Documents\ImmPanelRevis17\donnees\ipc_ocdeSML.dta", unm(m)
drop _merge
replace sn=sn*(100/ipc_ocde)

/* concatenate multiple jobs in one year */
save temp.dta, replace
contract nninouv an
gen mult=(_freq>1)
drop if _freq==1
rename _freq nmult
save mult.dta, replace
use temp.dta, clear
joinby nninouv an using mult.dta, unm(m)
drop _merge
drop if missing(mult)

/* fait la somme du nombre de jours 
et du salaire net sur les multiples emplois */
/* durée de la paye total */
/* durée paye normalisée & salaire réel total */
save tempCOL.dta, replace
collapse (sum) dps=dp sns=sn, by(nninouv an) fast
save dps.dta, replace
use tempCOL.dta, clear

/* fusionne avec */
/* durée de paye totale */
use temp, clear
joinby nninouv an using dps.dta, unm(m)
drop _merge

/* garde comme référence 
la firme (ie secteur) 
avec le plus grand salaire */
sort nninouv an sn
by nninouv an: drop if _n != _N
egen indiv=group(nninouv)
xtset indiv an
joinby nninouv an using mult.dta, unm(m)
drop _merge
replace mult=0 if missing(mult)

/* remplace salaire en salaire max et durée de
paie totale dans dp et sn pour ceux dont on a de multiples
activitées dans une année */
replace dp=dps if mult==1
replace sn=sns if mult==1
replace dp=360 if dp>360

gen fy=(dp==360)

/* Calcul du "revenu permanent" => moyennes autour des 5 RP */
/* utilise plusieur définitions */
save tempattrit.dta, replace

use tempattrit.dta, replace
egen time=group(an)
egen nindiv=group(nninouv)
drop if missing(nindiv)
keep if homme==1
gen first=(age>=25 & age<=50)
gen twice=(first==1 & !missing(F.time))
tab twice an if first, col
